Perturbation theory for confined systems 



Francisco M. Fernandez 
INIFTA (UNLP, CCT La Plata-CONICET), Division Quimica Teorica 
Diag. 113 y 64 (S/N), Sucursal 4, Casilla de Correo 16 
1900 La Plata, Argentina 
E-mail: fernande@quimica. unlp .edu. ar 

Abstract 

We discuss the application of perturbation theory to a system of 
particles confined in a spherical box. A simple argument shows that 
the particles behave almost independently in sufficiently strong con- 
finement. We choose the helium atom with a moving nucleus as a 
particular example and compare results of first order with those for 
the nucleus clamped at the center of the box. We provide a suitable 
explanation for some numerical results obtained recently by other au- 
thors. 

1 Introduction 

In a recent paper Montgomery Jr. et al |1] solved the Schrodinger equation 
for a He atom with its nucleus clamped at the origin of a box of radius Rc with 
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impenetrable walls. They applied perturbation theory for the case of strong 
confinement (sufficiently small Rc) and obtained the first five coefficients of 
the expansion (with different degrees of accuracy). One of the conclusions in 
that paper was that the interaction between the electrons decreases with the 
box radius. The authors illustrated this behavior by means of the overlap 
between the wavefunctions for the confined He"*" and for the free electron in 
the box. 

The purpose of this letter is to discuss those numerical results from a more 
general point of view. To this end in Sec. |2] we apply perturbation theory 
to a system of particles in a spherical box and discuss the behaviour of 
a more general overlap integral. As a particular example, we compare the 
energies (corrected through first order) of the He atom when the nucleus is 
clamped at the center of the box and when it moves confined in the box. 
Finally, in Sec. [3] we comment on the results and draw conclusions. 



2 Perturbation theory for strong confinement 

We consider a system of particles of masses rrii and charges g^. The 
nonrelativistic Hamiltonian operator is 

^2 TV V72 N-1 N 



where rjj = |rj — rj| is the distance between the pair of particles located at 
If the system is confined in a box of radius Rc with impenetrable walls, 



Ti and Tj 
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any solution ip to the time-independent Schrodinger equation 

H^p = Ell) (2) 

should vanish when > Rc for any given particle i. In order to apply per- 
turbation theory in the case of strong confinement -Rc — we first convert 
the Schrodinger equation ([2]) into a more convenient dimensionless eigenvalue 
equation. We choose a representative particle (say i = 1) and define dimen- 
sionless masses = mi/mi, charges = Qi/qi and coordinates r- = Vi/Rc 
(V- = i?cVj). We thus obtain a dimensionless Hamiltonian operator 

Rc Aneoh^ , . 

A = — , a = 2- (3) 

a rriiqf 

and the dimensionless eigenvalue equation 

Hdip = e^, e = -j^E = -^^E (4) 

The new boundary conditions are cp = if any r- > 1. Note that if the 
chosen reference particle is an electron, then mi = mg, qi = —e and a = ag 
is the Bohr radius. The transformation just proposed is a generalization of 
the one recently applied to the confined hydrogen atom [2]. 

It is clear that Hd{^ = 0) = is the dimensionless Hamiltonian operator 
for a system of free particles in a spherical box of unit radius. Therefore, 
we can solve the eigenvalue equation H^f^^^ = e^^^(f^^^ exactly in terms of 
products of spherical harmonics and Bessel functions [3]. It may also be 
necessary to consider the permutational symmetry of the wavefunction and 
add the corresponding spin factors [T]. 



For concreteness, let us consider the He atom. We assume that the parti- 
cles 1 and 2 are the electrons and the remaining one is the nucleus; that is to 



say: mi = 1712 = rrie and = rrin- Obviously, in such a case m'^ 



rrir, 



TTin/ rrie and the unperturbed wavefunction for the ground state is 



^sin(7rr^) sin(7rr2) sin(7rr3^ 
r[ r'2 rg 



Kl)/3(2)-/3(l)a(2)] (5) 



Note that the present model accounts for the motion of the nucleus and that 
the nuclear factor sin(7rr3)/r3 does not appear if this particle is clamped at 
the center of the box fHH]. 

If we apply straightforward Rayleigh-Schrodinger perturbation theory we 
obtain the well-known expansions 



j=o j=0 
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In particular, for the energy we have 
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that is a generalization of the result derived by Laughlin [5j and discussed 
by Laughlin and Chu [1] and Montgomery et al [1]. Note that equations ([6]) 
and ([7]) apply to the most general case of a system of N particles ([3]). 
If both Lf and if^^^ are normalized to unity we can easily prove that 



< 1, lim 



(0) 



(8) 



which clearly account for the behaviour of the overlap integral in Fig. 1 of 
Montgomery et al [Ij for the particular case of the He+. We stress that 
Eq. ([H]) applies to the general case of particles. 
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As an illustrative example we calculate the energy of the helium atom 
corrected through first order. When the nucleus is clamped at the origin we 
have (Hi] 

e(A) _ 9.8696044 7.9645404 

^ " A ^ ' 

On the other hand, when the nucleus moves the result is 

e(A) _ 9.870280744 5.358219501 

where we have chosen m'^ = 7296.300 mg. It is worth noting that the effect of 
the nuclear motion is more noticeable on the average potential energy than 
on the kinetic energy. [3]. 



3 Conclusions 

We have shown that converting the Schrodinger equation into a dimension- 
less eigenvalue one greatly facilitates the application of perturbation theory 
to strongly confined systems. In particular this approach clearly shows that 
the interaction between the particles becomes negligible as the confinement 
increases. In this way we could provide a suitable mathematical basis for 
recent numerical calculations on the He atom with a nucleus clamped at ori- 
gin. Eq. ([8]) not only explains the behavior of the overlap integral calculated 
by Montgomery et al [Ij but also reveals that the same kind of curve should 
be expected for any system of particles confined in a spherical box. 

We have also shown that the effect of the nuclear motion on the kinetic 
energy of the confined atom is not as important as its effect on the average 
of the Coulomb interactions. The reason is that the energy of the confined 
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atom changes markedly with the location of the clamped nucleus. Therefore, 
when it moves there is a sort of average contribution to the potential energy 
from all the possible locations inside the box. 
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